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Abstract 

The  ultrasonic  plane  wave  reflection  coefficient  for  a  flat  solid  plate  immersed  in  a  liquid 
should  be  reducible  to  that  for  a  liquid- solid  interface  as  the  plate  thickness  becomes  very 
large,  hxisting  formulae,  which  are  based  on  the  assumption  that  propagation  vectors  are 
real  quantities,  do  not  provide  this  reducibility.  It  is  shown  that  reducibility  can  be  achieved 
only  when  the  propagation  vectors  are  complex  quantities,  i.e.,  when  absorption  is  taken 
into  account.  The  approach  given  here  can  thus  be  used  to  determine  reflectivity  for 
nutterial  thicknesses  which  are  neither  truly  infinite  nor  very  thin.  The  results  provide 
practical  criteria  to  determine  the  lower  limit  of  sample  thickness  that  can  be  considered 
a  half  space. 


Introduction 

Reflection  of  ultrasonic  plane  waves  in  a  liquid- solid-liquid  system  has  been  described  1 1 .2) 
by  the  amplitude  reflection  coefficient  RLSL .  derivable  from  consideration  of  the  bound¬ 
ary  conditions  at  the  interfaces.  The  plane  wave  amplitude  reflection  coefficient  RLS  for 
a  liquid-solid  interface  was  also  obtained  (1.2)  when  the  boundary  conditions  at  that 
interface  are  taken  into  account. 

It  was  shown  (2)  that,  beside  the  material  characteristics,  the  behavior  of  RLSL  is  solely 
determined  by  one  physical  parameter,  the  product  frequency  f  times  the  thickness  d  of 
the  solid  plate.  Mathematical  and  physical  consistency  requires  that  Rlsl  become  RLs  as 
this  product  tends  to  infinity.  This  property  has  not  been  established,  although  it  was 
briefly  discussed  by  Pitts  [2J.  It  is  also  interesting  to  understand  the  behavior  of  RLsl 
in  this  limit  as  well  as  to  be  able  to  determine  the  thickness  of  a  solid  plate  that  can  be 
considered  a  practically  acceptable  half-space. 

Preliminary  investigation  indicates  that  absorption  plays  an  essential  part  when  Rlsl 
examined  in  the  limit:  frequency  x  thickness  The  present  paper  will  study  the  exist¬ 
ence  of  such  a  limit  with  and  without  absorption  being  taken  into  consideration.  As  a 
result,  reflectivity  characteristics  of  a  thick  solid  plate  will  also  be  provided. 


1.  Plane  wave  reflection  coefficients 

The  reflection  coefficient  Rlsl  for  a  solid  immersed  in  a  liquid,  and  that  for  a  liquid- 


solid  interface.  Rls-  are  given  by  |2] 

Rls  =  (N,  +N, -N,)/(N,  +N2  +N3).  (1) 

Rlsl  =  N  (FaFs).  (2) 

where 

N,  =  <k,J  -  2kJ  t:.  (21 

N2  =  4kjKjK».  (4) 

N3  -  pkjKd/K,  (5) 

N  *  Nj  +  N$  -  Nj  +  2N|  N2  ( 1  -  cosPcosQ)  /  (sinPsinQ).  (61 
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Ft  =  |(  1  -  cosP)  /  sinPJ  N'i  +  |(  1  -  cosQ)  I  smOJ  Nj  +  iN3.  (7) 

Fs  =  |(  1  +  cosP)  /  sinP]  N,  +  |(  1  +  cosQ)  /  sinQ|  N2  -  iN3.  (8) 

The  symbols  used  in  Eqs.  ( 3-8)  are  defined  as  follows: 

kx  =  ksinOj.  (y) 

k  =  2trf/v 

kd  =  2trf/vd  without  absorption.  (10) 

ks  =  2»rf/v5 

k  =  (2fff/v)(  1  +  ia/2»r) 

kd  =  (2trf/vd)(l  +  iad/2rr)  with  absorption  (II) 

ks  =  (2rrf/vs)  ( I  +  ias/2rr) 

K  =  (ks-kj Kd  =  (k^-k’)*.  Ks  =  (kf- kj (12) 

P  =  dKd,  0  =  dKs,  (13) 

where 

=  angle  of  incidence, 
p  =  liquid  density  /  solid  density. 


and  where  v,  vd.  vs  and  a.  ad,  a,  are  the  measured  sound  velocities  and  attenuation  per 
wavelength  of  the  longitudinal  waves  in  the  liquid,  and  the  longitudinal  and  shear  waves 
in  the  solid,  respectively.  Consideration  of  Eqs.  (1-13)  shows  that  Rlsl  depends  on  Id 
only  through  terms  involving  P  and  Q.  The  existence  of  lim  Rlsl  thus  dwells  essentially 
on  the  behavoiur  of  P  and  Q  at  such  a  limit.  fd*” 


II.  Limit  of  Rlsl  as  fd  tends  to  infinity 

When  absorption  is  included,  that  is  a,  ad ,  a5  =£  0,  P  and  Q  have  the  lorm 
P  =  fd(pR  +  ipi)  and  Q  =  fd(qR  +  iq( ). 

where  pR,  pt,  qR ,  and  q(<  in  general  are  non-zero  and  finite.  The  form  of  the  wave  func¬ 
tions  chosen  in  the  derivation  of  Rlsl  implies  |2]  that  pi  and  qi  are  non-negative,  so  as 
to  make  the  sound  amplitude  decay  as  physically  expected.  In  the  limit  of  fd  tending  to 
infinity,  the  behaviour  of  RLsl >s  governed  by  terms  containing  P  and  Q.  It  is  observed 
that 


cosP  =  cos(fdpR  )  cosh(fdpr )  -  isin(fdp«  )  sinh(fdpi ). 
sinP  =  sin(fdpR  )  cosh(fdpj )  +  icos(fdpR  )  sinh(fdpi ). 
Since  p(  is  positive 

lim  cosh(fdpi)  =  lim  sinh(fdpi)  =  lim  ^e^P*. 
fd-***  fd-***  fd-~  - 

Hence 


(14) 


(15) 
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Hence 


lim  cosP  =  lim  i  e  dPi  [cos(fdpR  )-  isin(fdpR)], 
f'4  fd-»  *■ 


flim  sinP  =  ^im  ^ie^P*  |cos(fdpR  )  -  isin(fdpK )]  • 


(16) 


bquations  ( 14-16)  are  now  used  to  examine  the  limit  when  fd-»°°  of  terms  involving  P 
and  0  in  the  expression  of  RLSL  First  consider 


l±^efdPi  |cos(fdpK  )  -  isin(fdpR  )] 
fd»-  sinP  i,efdPi  Jcos(fdpK  )  —  isin (fdpK  )J 


lim 


itcosP 


=  lim 


but 


fhm  1 1  +  ie^P'cosddpn  )J  =  fhm  je^Pi  [e*^dPi  +  iCOs(fdpR  ))} . 

=  =  ^Ijm  {±ie^dP*cos(fdpK )J , 

since  p!  is  always  positive,  as  mentioned  earlier. 

Hence 

<l±cosP)  ±ieldP<  |cos(fdpR)-isin(fdpR)J 

Inn  - - —  =  lint  — -  . 

rd-  s.nl  4iefdPi  |Cos(fdpR)-isin(fdpR)] 


(17) 


(IK) 


(19) 


lim 

fd»“ 


(  1  tcosP) 
sinP 


Similarly,  it  can  be  shown  that 


+  i. 


lim 

fd-- 


( 1  ±  cosQ) 
sinQ 


+  i. 


Next  consider  the  same  limit  of  the  remaining  term  involving  P  and  0, 


lim  1-cosPcosQ 
sinPsinQ 


In  view  ol  bq.  ( 16).  the  real  part  of  ( 1-cosPcosQ)  is  seen  to  behave  as 


(20) 


(21) 


lim  Re ( 1-cosPcosQ)  = 

W  - 

lim  ||_Iefd(Pi +tli)  (Cos(fdpR  )cos(fdqK  )-sin(fdpR  )sin(fdqR  ))} , 
fd-  4  '  (22) 

=  la™  j  -  ^e*^d^Pi  +9i  1  (cos(fdpR  )cos(fdqR  )-sin(fdpR  )sin(fdqR )]} , 
where  a  manipulation  similar  to  that  of  Eq.  (18)  has  been  performed. 


Use  of  Eq.  (22)  yields 

..  1-cosPcosQ 
lim  — — — ~ 
fd--  smPsmQ 

-le^Pt  +  9tl(cosfdpR  -  isin  fdpK)(cosfdqR- isin  fdqR) 

=  lim  j-T  — - - -  . 

w*~ i  i3 efd(P»  +  91 )  (cosfdpR  -isinfdpR  )  (cosfdqR  -  i sin fdqH ) 

or 

.  1-cosPcosQ 

lim  - - — — -  =  1. 

fd»«  sinPsinQ 


(23) 


(24) 
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Substitution  of  Eqs.  (20), (2 1)  and  (24)  into  Eq  (2)  results  in 

_  N?+Nj-N|+2N,N, 

fd™  lsl  ‘  (jNi  +  iNj  +jN3)(_iN)  _iNj  _jN3) 

_  (N,  +  Nj )5  -  N|  =  N,  +Nj  -N3 
"  (N,  +  Nj  +N3)5  N,  +N5  +N, 


or 


lim  Rlsl  =  lim  Rls. 

fd  -  •»  fd-»- 


(25) 

(26) 

(27) 


III.  Discussion 

The  above  derivation  establishes  that  the  plane  wave  reflectivity  from  a  solid  plate  im¬ 
mersed  in  a  liquid  will  become  that  from  a  liquid-solid  interface  when  the  plate  thickness, 
normalized  with  respect  to  the  acoustic  wavelength,  tends  to  infinity.  This  result  holds 
only  when  absorption  in  the  media  is  taken  into  account.  Without  absorption  being  con¬ 
sidered,  such  limit  does  not  exist.  In  particular,  when  absorption  is  neglected  it  can  be 
shown  that  RLsl  tends  to  RLS  only  when  the  incident  angle  is  larger  than  the  shear  crit¬ 
ical  angle.  This  is  graphically  illustrated  in  Figure  1 ,  where  the  modulus  of  Rlsl  is  plotted 
against  the  incident  angle  with  a  very  high  value  of  fd,  50km/s,  for  a  water- brass- water  sys¬ 
tem.  There  the  IRlsl  '  curve  is  seen  to  be  irreducible  to  the  IRLs  I  °ne  when  absorption  is 
neglected.  It  is  worth  noting  that  near  the  Rayleigh  critical  angle  Rlsl  always  tends  to 
RLs  n«  matter  whether  absorption  is  considered  or  not.  This  confirms  the  understanding 
that  Rayleigh  surface  waves  or  Rayleigh-type  Lamb  waves  are  always  physically  confined 
to  the  first  interface,  which  is  the  one  on  which  the  incoming  sound  waves  are  incident. 
The  physical  location  of  the  second  interface  obviously  does  not  affect  the  behaviour  of 
these  propagation  modes. 


10  20  30  40  SO 

ANGLE  OF  INCIDENCE  C) 


Figure  1.  Reflection  amplitude  for  a  lossless  water-brass-water  system  with  fd  '  50km/s. 
The  solid  curve  is  reflection  amplitude  for  a  lossless  water-brass  system. 
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ire  2.  Reflection  amplitude  for  an  absorptive  water-brass-water  system  with  fd  = 
31.75km/s  and  f  *  5MHz.  The  reflection  amplitude  for  a  water-brass  inter¬ 
face  with  absorption  for  the  same  frequency  is  illustrated  by  the  solid  curve. 


Figure  3.  Reflection  amplitude  for  an  absorptive  water-brass-water  system  with  f  = 
10MHz  and  (a)  fd  =  31.75km/s  and  (b)  fd  =  63.5km/s,  respectively.  The  reflec¬ 
tion  amplitude  for  absorptive  water-brass  systems  of  the  same  frequency  is 
plotted  in  the  solid  curve. 
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Figures  2  and  3  demonstrate  graphically  the  behaviour  of  IRLSL  I  when  absorption  is 
included  for  a  water-brass- water  system.  Attenuation  data  used  for  brass  are  derived  from 
measurements  made  by  Papadakis  |3]  and  those  for  water  provided  by  Pinkerton  |4] . 
The  modulus  of  Rlsl  is  presented  in  Figure  2  for  f  =  5MHz.  fd  =  3 1 .75km/s.  and  in 
Figure  3  for  f=  10MHz.  fd  =  31 .75  and  63.5km/s.  Comparison  of  the  three  curves  of 
Figures  2  and  3  leads  to  the  following  observations:  (i)  as  fd  becomes  larger  IRLsl  I 
behaves  more  like  IRLs  k  and  (ii)  for  the  same  fd  value  IRlsl  I  tends  to  IRls  I  faster  for 
higher  absorption. 


Figure  4.  Reflection  amplitude  for  an  absorptive  water-Plexiglas-water  system  with  (a) 
fd  =  3.18km/s,  f  *  2MHz  and  (b)  fd  =  1 2.7km/s,  f  *  4MHz.  The  reflection  amp¬ 
litude  for  an  absorptive  water-Plexiglas  interface  at  2MHz  is  described  by  the 
solid  curve. 


Figure  5.  Reflection  amplitude  for  an  absorptive  water-Plexiglas-water  system  with  fd  = 
6.35km/s  but  (a)  f  =  2MHz  and  (b)  f  =  4MHz.  The  reflection  amplitude  for  an 
absorptive  water-Plexiglas  interface  is  described  by  the  solid  curve. 


Since  the  physical  properties  of  Plexiglas  are  remarkably  different  from  brass,  it  is  int¬ 
eresting  to  see  how  IRlsl  I  behaves  in  comparison  with  IRls  !  for  a  water-Plexiglas-water 
system.  Similar  results  are  obtained  in  this  case  except  for  the  fact  that  with  a  highly 
absorptive  substance  such  as  Plexiglas  RLsl  reduces  to  Rls  at  much  lower  values  of  fd. 
These  are  illustrated  in  Figures  4  and  5.  The  attenuation  data  used  for  Plexiglas  are  taken 
from  Ngoc  (S) .  The  same  observations  as  made  for  water-brass-water  also  apply  here: 
IRlsl  I  approaches  IRls  I  faster  for  higher  fd  values  and  absorption. 
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IV.  Conclusion 

The  above  results  provide  a  unified  picture  of  plane  wave  reflectivity  for  liquid-solid  and 
liquid-solid-liquid  layered  systems.  It  is  seen  that  absorption  plays  an  essential  role  in 
verifying  that  RLSL  >s  reducible  to  RLS  when  the  thickness  of  the  solid  plate  becomes 
very  large  compared  to  the  acoustic  wavelength.  The  question  of  how  thick  a  solid  plate 
can  be  considered  to  be  a  half-space  for  practical  purposes  must  be  addressed  with  proper 
consideration  of  the  absorption  characteristics  of  the  media  involved. 


Acknowledge  mem 

Part  of  this  work  was  supported  by  the  Office  of  Naval  Research.  U.S.  Navy.  Two  of  the  authors 
(O.J.L.  and  W.G.M.)  gratefully  acknowledge  support  through  a  Grant  of  the  Scientific  Affairs  Division 
of  NATO,  which  made  cooperation  between  the  authors  possible. 


References 

(I]  Brekhovskikh,  L.M.,  Waves  in  Layered  Media,  Academic  Press,  New  York,  I960. 

[2J  Pitts,  L.E.,  A  Unified  Theoretical  Description  of  Ultrasonic  Beam  Reflections  fora 
Solid  Plate  in  a  Liquid,  Ph.D.  Thesis,  Physics  Department,  Georgetown  University, 
Washington,  D.C.  1975. 

[3]  Papadakis,  E.P.,  J.  Acoust.  Soc.  Am.,  37, 71 1,  1965. 

|4]  Pinkerton, J.M.M.,  Nature,  160,128.1947. 

15  J  Ngoc,  TDK.,  Influence  of  Absorption  on  Ultrasonic  Nonspecular  Reflectivity, 

Ph.D.  Thesis,  Physics  Department,  Georgetown  University,  Washington,  D.C. ,  1979. 

(Received  1 9  October  1981) 


Acoustics  Letters  VoL  S,  No.  3,  1981 


112 


IEEE  TRANSACTIONS  ON  SONICS  AND  ULTRASONICS,  VOL.  SU-29,  NO.  2,  MARCH  1982 


Correspondence 


Influence  of  Plate  Mode  Structure  and  Gaussian  Beam  Profile 

Characteristics  on  Ultrasonic  Reflection  and  Transmission 

TRAN  D.  K.  NGOC  and  WALTER  G.  MAYER 

AJa tract -It  is  shown  why  ultrasonic  beam  reflection  from,  and 
transmiasion  through,  a  solid  plate  immersed  in  a  liquid  may  be  non- 
apecular,  depending  cm  beam  width  and  the  structure  of  the  normal 
modes  of  vibration  of  the  plate.  The  analysis  is  carried  out  for  Gaus¬ 
sian  incident  beams. 

INTRODUCTION 

The  profile  of  an  acoustic  beam  reflected  from  or  trans¬ 
mitted  through  a  layered  medium  has  recently  been  calculated 
by  several  authors.  An  analytical  method  was  used  by  Bertori 
and  Tamir  [  1  ]  for  a  liquid-solid  interface  and  by  Pitts  [  2]  and 
Ng  [31  for  a  solid  plate  immersed  in  a  liquid,  respectively. 

Ngoc  and  Mayer  (4]-[61  determined  these  profiles  numeri¬ 
cally  and  showed  that  the  use  of  a  numerical  integration 
method  eliminated  some  of  the  limitations  encountered  in  the 
analytical  formulation. 

For  a  solid  plate  immersed  in  a  liquid,  Ngoc  and  Mayer  [6] 
investigated  the  general  features  of  the  reflected  and  trans¬ 
mitted  beam  profiles  for  incidence  at  plate-mode  and  between- 
mode  angles.  Hate-mode  angles  of  incidence  are  defined  as  the 
incident  angles  of  the  impinging  plane  waves  at  which  leaky 
Lamb  waves  are  generated,  while  the  latter  type  of  incidence 
denotes  an  angle  of  incidence  located  exactly  between  two 
adjacent  plate-mode  angles.  Changes  in  profiles  were  also 
studied  as  a  function  of  angle  of  incidence  and  beamwidth. 
However,  there  are  neither  systematic  studies  of  profile 
changes  nor  of  the  relationship  between  the  mode  structure  of 
a  loaded  solid  plate  and  the  beam  structure  of  a  Gaussian 
beam. 

This  correspondence  is  concerned  with  profile  changes  in  a 
broader  perspective  in  order  to  provide  better  practical  insight 
into  nonspecular  reflection  and  transmission  phenomena.  A 
nonspeculariy  reflected  or  transmitted  beam  will  be  viewed  as 
a  result  of  the  interplay  between  plate-mode  structure  and 
beam  characteristics.  Emphasis  will  therefore  be  placed  on 
using  an  adequate  range  of  values  for  beamwidth  as  well  as  a 
more  representative  mode  structure. 

Mathematical  and  Physical  Description 

The  sound  field  of  a  bounded  beam  reflected  from  or  trans¬ 
mitted  through  a  solid  plate  immersed  in  a  liquid  is  repre¬ 
sented  (6)  by 

Uw(x,  0)  *  (27t)-‘  1  R(kx)  V(kx) 

*'*/-»  /w, 

•exp(f**x)d*„  (1) 

rkt**iw, 

f/(7V,  -d)  -  (2ir)-‘  I  T(/tx)  V{kx ) 

•expli(*,x-M)l<**x.  (2) 
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where  U ^  and  are  the  sound  amplitudes  of  the  re¬ 
flected  and  transmitted  fields,  and  R(kx)  and  T(kx)  are  the 
plane  wave  reflection  and  transmission  coefficients.  The  x- 
direction  is  given  by  the  propagation  direction  of  the  Lamb 
wave,  and  the  thickness  of  the  plate,  d,  is  measured  along  the 
z-direction.  The  other  symbols  used  in  above  expressions  are 
defined  as  follows: 

kx  x -component  of  the  wave  vector  k, 

k ,  =  (fc2-k2)>/2, 

V(kx)  Fourier  transform  of  incident  beam,  defined  in  (4), 
kt  =  (uijv)  sin  0(, 

to  sound  angular  frequency, 

6/  incident  angle, 

v  sound  speed  in  the  liquid, 

W0  =  W/cos  6it 

W  half-width  of  the  Gaussian  beam. 

Equations  (1)  and  (2)  are  based  on  the  fact  t  n  incident 
beam  can  be  represented  by  a  superposition  ol  .•  ite 
number  of  plane  waves  of  different  incident  an  '  ,ne  re¬ 
flected  and  transmitted  beams  are  then  constructed  from 
the  constituent  plane  waves  after  they  have  interacted  with 
the  layered  system  via  the  reflection  or  transmission  coeffi¬ 
cient,  derived  from  the  plane  wave  theory.  In  general,  the  re¬ 
sulting  reflected  or  transmitted  beam  profiles  are  determined 
by  the  physical  characteristics  of  the  incident  beam  and  the 
layered  medium,  which  here  is  taken  to  be  a  solid  plate  im¬ 
mersed  in  a  liquid. 

A.  Beam  Characteristics 

A  Gaussian  beam  of  width  2W  incident  at  an  angle  0,-  is 
characterized  at  z  =  0,  the  boundary  between  the  liquid  and 
the  solid  plate,  by 

tfinc(*.  0)  =  «P  l~(x/W0)2  +  ixk/),  (3) 

whose  Fourier  transform  in  the  kx  domain  is 

V(kx)  =  i r>/»  w0  exp  [-(**  -  *,)’(FtW2)2 ].  (4) 

Equation  (3)  implies  that  the  incident  angles  of  constituent 
incoming  plane  waves  are  perturbed  about  0/  and  that  this 
angular  spreading  is  a  function  of  the  beamwidth.  Fot  a  given 
sound  frequency,  this  angular  spreading  can  be  measured  in 
terms  of  0^,  an  angle  in  degrees,  expressing  the  half-width  of 
the  Gaussian  function  V(kx ),  and  0hw  can  be  written  [6]  as 

0HW*sm-‘  (sin  0{  +  0.53  Wk)  -  sin*1  (sin  0,  -  0.53  W*), 


where  Wk  -  ( nv  cos  0/)/(wRO. 

Some  simple  manipulation  of  (5)  gives  rise  to 

0HW-3O.4X/H',  (6) 

where  X  is  the  sound  wavelength  in  the  liquid.  The  identity 
(6)  is  a  good  approximation  when  \/W  is  less  than  unity. 
Typical  values  of  0hw  416  summarized  in  Table  I  for  a  fre¬ 
quency  of  2  MHz  and  for  several  values  of  W  under  considera¬ 
tion. 

B.  Mode  Structure 

The  transmission  or  reflection  coefficient  for  a  solid  plate 
immersed  in  a  liquid  can  be  used  to  describe  the  acoustical 
characteristics  of  the  plate.  In  general,  these  coefficients  are 
functions  of  the  product  fd  (frequency  times  plate  thickness). 
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TABLE  l 

Angular  Spreading  of  a  Gaussian  Incident  Beam  for 
Various  Beamwidth  and  Frequency  of  2  MHz 


w 

(in) 

K/W 

®HW 

(degrees) 

i 

0.2362 

14.4 

i 

0.1181 

7.2 

i 

0.0591 

3.6 

l 

0.0394 

2.4 

2 

0.0147 

0.89 

8 

0.0037 

0.22 

32 

0.0009 

0.06 

TABLE  II 

Mode  Structure  of  a  Water-Brass-Water  System  for 
fd  =  7  MHz  mm 


eP 

Modes  (degrees)  S(p„ ,  pm ) 


-40(So> 

43.2 

39.6 

3.6 

S , 

34.2 

5.4 

M 

27.4 

6.8 

s2 

22.05 

5.35 

'U 

19.3 

2.75 

s3 

16.25 

3.05 

s4 

16.0 

0.25 

a4 

15.25 

0.75 

Ss 

5.65 

9.6 

S,  4S,A,S,  a,  s,  A,  A. 


A  !  ! 


I 


!  I  I  i 


:  5 


0  10  20  30  AO  50 

INCIDENT  ANGLE 


Fig.  1.  Plane  wave  transmission  modulus  for  water-brass-water  system 
with  fd  =  7.0  MHz  •  mm.  Standard  plate-mode  nomenclature  is  used 
to  identify  plate-mode  angles  of  incidence. 


Their  analytical  behavior  is  completely  determined  by  a  num¬ 
ber  of  pole-zero  pairs  [6].  These  pole-zero  pairs  are  associated 
with  plate-mode  angles  of  incidence.  The  complete  set  of 
plate  modes  forms  the  mode  structure  of  the  layered  system 
under  consideration.  Thus  the  mode  structure  identifies  the 
plate-mode  angles  dp  and  the  angular  separation  between 
them.  In  the  following,  S{p„ ,  pm  )  denotes  the  angular  separa¬ 
tion  between  modes p„  and  pm.  Table  II  describes  the  mode 
structure  of  a  brass  plate  immersed  in  water  for  an  fd  of  7 
MHz  -  mm.  This  value  of  fd  was  chosen  so  that  S(pn,  pm )  has 
a  wide  range  of  values,  which  are  also  comparable  to  the  values 
of  0hw  listed  in  Table  I.  The  corresponding  transmission  co¬ 
efficient  is  shown  in  Fig.  1.  In  both  Table  II  and  Fig.  1,  the 
letters  A  and  S  denote  symmetric  and  antisymmetric  modes, 
respectively. 


Calculated  Results 

The  exptessions  given  above  can  be  used  to  calculate  the 
energy  or  amplitude  profiles  of  reflected  or  transmitted  sound 
fields  produced  by  an  incident  Gaussian  beam.  It  is  well- 
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Fig.  2.  Transmitted  and  reflet,  led  beam  profiles  for  water-brass-water, 
7  MHz  •  mm  and  plate-mode  incidence  at  39.6°.  Beamwidth  values 
are  (a)  |  in,  (b)'J  in,  (c)  2  in,  and  (d)  8  in. 


known  that  the  reflected  or  transmitted  sound  field  quite  fre¬ 
quently  is  no  longer  Gaussian,  that  the  reflected  beam  often  is 
“displaced”  [7],  sometimes  split  into  distinct  beams  with  an 
energy -null  separating  the  sections  [8],  or  sometimes  reflected 
as  a  Gaussian  beam  without  any  apparent  change  from  the 
original  Gaussian  profile. 

These  nonspecular  reflection  and  transmission  phenomena 
can  be  predicted  if  one  takes  into  account  the  mode  structure 
of  the  plate,  the  width  and  frequency  of  the  incident  beam, 
and  the  angle  of  incidence  as  given  in  (l)-(6). 

As  an  example  of  sue  calculations,  consider  a  brass  plate 
immersed  in  water,  with  file  parameter  fd  bein’  7  MHz  -  mm 
and  the  Gaussian  incident  beam  having  a  frequency  /  =  2  MHz. 
Two  series  of  calculations  are  discussed  here  for  two  types  of 
incidence:  incidence  at  the  plate-mode  angle  0,  =  39.6°  and 
incidence  at  the  u? tween-mode  angle  0,-  =  4 1 .4°. 

For  the  selected  A  i  plate-mode  incidence  angle,  four  pro¬ 
files  are  calculated  for  the  following  beamwidths:  (a)  |  in, 

(b)  ^  in,  (c)  2  in,  and  (d)  8  in.  The  resulting  profiles  are 
plotted  in  Fig.  2. 

(a)  Since  |  in  is  a  rather  narrow  beamwidth,  the  angular 
spreading  0hw  covers  .ree  adjacent  modes,  A0  on  the  left 
and  S |  and  A2  on  the  .jht.  With  this  wide  coverage  of  0hw. 
the  reflected  profile  sh  ws  a  complicated  pattern  with  more 
than  one  secondary  peaks  and  intensity  nulls,  caused  by  the 
influence  of  the  involved  neighboring  modes.  The  transmitted 
beam  exhibits  a  lateral  displacement  of  approximately  2 WQ, 
which  is  the  projection  of  the  beamwidth  onto  the  reflecting 
interface. 

(b)  A  beamwidth  of  ^  in  implies  that  the  angular  spreading 
0HW  is  7.2°,  including  on  both  sides  two  modes,  A0  and  S j . 
The  reflected  profile  exhibits  only  one  intensity  null  and  the 
only  secondary  peak,  broadened  considerably,  indicates  com¬ 
peting  influences  of  the  adjacent  modes.  In  the  transmitted 
profile,  this  phenomenon  also  appears  in  the  trailing  edge  of 
the  sound  field,  while  the  lateral  displacement  is  decreased  but 
still  significant,  extending  to  approximately  W0. 

(c)  For  a  beamwidth  of  2  in,  the  angular  spreading  0hw  = 
0.89°  is  narrowly  centered  about  the  A  j  mode.  Both  reflected 
and  transmitted  profiles  are  influenced  only  by  this  mode. 
These  profiles  are  typical  whenever  only  one  mode  is  involved. 

(d)  With  a  very  large  beamwidth,  8  in  in  this  case,  the  angu¬ 
lar  spreading  becomes  extremely  small  such  that  the  beam 
profiles  are  strongly  determined  by  the  behavior  of  the  plane 
wave  coefficients.  Almost  total  transmission  is  obtained  with 
no  nonspecular  features  present.  The  reflected  profile,  al¬ 
though  of  very  low  intensity,  retains  its  nonspecular  features 
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Fig.  3.  Transmitted  and  reflected  beam  profiles  for  water-brass-watei. 
7  MHz  ■  mm  and  between-mode  incidence  at  41.4°.  Beamwidth  value 
values  are  (a)  $  in,  (b)  2  in,  (c)  i  in,  and  (d)  $  in. 


because  of  an  extremely  sharp  phase  shift  of  the  reflection  co¬ 
efficient  at  plate-mode  incidence. 

The  results  of  calculations  for  the  between-mode  angle  of 
incidence  are  shown  in  Fig.  3  for  beamwidths  of  (a)  i  in, 

(b)  2  in,  (c)  ^  in,  and  (d)  ^  in.  The  following  features  should 
be  noted. 

(a)  The  beamwidth  of  ^  in  was  chosen  as  the  first  calcula¬ 
tion  value  for  a  specific  reason.  Special  nonspecular  features, 
consisting  of  a  series  of  closely  spaced  secondary  peaks,  were 
observed  to  be  most  pronounced  for  this  beamwidth  in  both 
reflected  and  transmitted  profiles.  In  this  case  0HW  =  2.4°, 
which  is  of  the  right  magnitude  to  induce  competing  influ¬ 
ences  of  modes  A0  and  A  ( . 

(b)  When  the  beamwidth  is  much  larger,  2  in,  then  0HW  = 
0.89°.  This  is  small  enough  to  eliminate  essentially  the  in¬ 
fluence  of  any  mode;  thus  nonspecular  phenomena  disappear 
and  both  the  high  reflection  and  low  transmission  are  entirely 
specular. 

(c)  For  a  beamwidth  of  ^  in,  the  angular  spreading  becomes 
rather  large,  with  0HW  =  7.2°.  This  includes  only  the  modes 
A0  and  A  |,  and  the  resulting  transmitted  and  reflected  pro¬ 
files  indicate  that  a  combination  of  two  modes  behaves  like  a 
single  mode  if  the  value  of  0hw  of  the  incoming  beam  is 
much  larger  than  the  mode  spacing. 

(d)  With  the  beamwidth  being  |  in  and  an  associated  0HW 
of  3.6°,  the  resulting  profiles  are  a  mixture  of  those  described 
in  (a)  and  (c)  above,  manifesting  both  interference  and  single¬ 
mode  nonspecular  features  in  the  reflected  and  transmitted 
profiles. 

Conclusion 

It  has  been  demonstrated  that  the  basic  features  of  non- 
specularly  reflected  and  transmitted  beam  profiles  for  a  solid 
plate  immersed  in  a  liquid  can  be  predicted  from  comparison 
of  the  incident  beamwidth,  i.e.,  its  angular  spreading  of  its 
Fourier  transform  and  the  mode  structure  in  the  neighborhood 


of  the  incident  angle.  The  observations  can  be  summarized  as 
follows. 

1)  Plate-Mode  Incidence:  For  a  rather  wide  Gaussian  beam, 
i.e.,  very  small  angular  spreading,  nonspecular  features  are  not 
significant,  and  beam  profiles  are  determined  essentially  by 
plane  wave  reflection  coefficients.  When  the  angular  spreading 
covers  several  modes,  i.e.,  for  narrow  beams,  the  beam  profiles 
become  nonspecular  and  complex,  indicating  the  competing 
influence  of  these  modes.  Typical  single-mode  profiles  are 
obtained  when  the  angular  spreading  of  the  beam  is  significant 
and  includes  only  one  mode. 

2)  Between-Mode  Incidence:  Special  nonspecular  phenom¬ 
ena  are  observable  when  0HW  —  S(p„,  pm)!2.  When  the  angu¬ 
lar  spreading  is  very  small,  reflected  and  transmitted  beam 
profiles  are  specular.  For  a  very  large  angular  spreading,  the 
influence  of  two  modes  is  essentially  the  same  as  that  of  a 
single  mode,  and  the  resulting  beam  profiles  are  similar  to 
single  mode  profiles. 

The  analysis  presented  here  can  also  be  used  to  explain  the 
results  of  Firestone  ( 9],  who  was  unable  to  produce  total 
transmission  through  plates  at  all  the  angles  of  incidence  where 
a  plate  mode  should  have  been  generated;  neither  S(np,  pm ) 
nor  0hw  were  considered  in  his  experiments. 

In  view  of  the  approach  outlined  above,  one  can  also  explain 
why  Schoch  |7|  was  unable  to  identify  experimentally  all 
possible  modes  of  a  plate  of  a  given  fd.  His  criterion  for  the 
excitation  of  a  mode  was  the  observation  of  “lateral  beam 
displacement”;  however,  such  a  displacement  occurs  only 
when  the  beamwidth  and  the  mode  spacing  have  a  very  spe¬ 
cific  relationship  to  each  other.  Using  Schoch’s  criterion, 
modes  other  than  those  identified  by  him  can  be  observed 
simply  by  adjusting  the  beamwidth  in  relation  to  the  known 
mode  spacing. 

Finally,  applying  the  present  analysis  to  the  known  mode 
spacing  where  fd  becomes  very  large  (the  typical  liquid-solid 
single  boundary )  one  can  see  why  only  one  of  the  extremely 
large  number  of  possible  modes,  the  Rayleigh  mode,  will  yield 
nonspecular  reflection.  As  fd  becomes  very  large,  all  the  pos¬ 
sible  Lamb  modes  are  extremely  closely  spaced  and  only  the 
lowest  mode,  the  Rayleigh  mode,  is  angularly  separated  from 
the  rest.  Thus  even  a  rather  wide  incident  beam  has  a  value  of 
0HW>  which  although  it  may  well  be  much  less  than  1°,  is 
much  larger  than  the  S(pn,  pm  )  for  the  Lamb  modes,  and  as 
was  pointed  out  above,  this  situation  will  result  in  simple 
specular  reflection.  Only  the  separated  Rayleigh  mode  is  able 
to  produce  nonspecular  reflectivity  effects,  provided  the  inci¬ 
dent  beam  has  the  correct  angular  spreading. 
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Characterization  of  Localized  Surface  Elastic  Defects 
by  Nonspecular  Reflection  at  the  Rayleigh  Angle 
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The  relationships  between  characteristics  of  elastic  defects  and  nonspecular  features  of 
bounded  ultrasonic  beams  reflected  at  the  Rayleigh  angle  from  a  liquid-solid  interface  are 
investigated.  The  results  can  serve  as  a  theoretical  basis  for  interpretation  of  Rayleigh  angle 
nonspecularly  reflected  beam  profiles  as  characterization  of  localized  surface  elastic  defects. 


KEY  WORDS:  bounded  ultrasonic  beam  reflectivity;  Rayleigh  angle:  beam  profile;  surface  defects; 
nonspecular  reflection:  NDE. 


1.  INTRODUCTION 

When  an  ultrasonic  bounded  beam  is  incident  at 
a  liquid-solid  interface,  the  reflected  beam  very  often 
has  a  lower  intensity  level  but  the  same  intensity 
distribution  profile.  However,  at  Rayleigh-angle  inci¬ 
dence,  the  reflected  intensity  profile  is  known  to 
display  a  number  of  distinct  features  that  substan¬ 
tially  deviate  from  the  profile  of  the  incident  beam. 
These  features,  which  are  normally  referred  to  as 
nonspecular  reflection  effects,  may  include  a  lateral 
beam  displacement,  one  or  more  minimum  intensity 
areas,  and  a  trailing  sound  field.  Such  a  nonspecu¬ 
larly  reflected  profile  is  illustrated  in  Fig.  1. 

Several  authors/1-2'  although  using  different 
mathematical  techniques,  were  able  to  describe  non¬ 
specular  reflectivity  at  the  critical  Rayleigh  angle  for 
a  liquid-solid  interface.  Nonspecular  reflectivity  was 
also  demonstrated  and  described*3'  for  incidence  near 
the  longitudinal  critical  angle.  Study  of  nonspecular 
effects  was  then  extended  to  a  more  complex  struc¬ 
ture,  that  of  a  solid  plate  immersed  in  a  liquid,  in 
which  both  reflected  and  transmitted  beams  were 
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Fig.  t.  Coordinate  system  and  a  typical  nonspecularly  reflected 
beam  profile  at  the  Rayleigh  angle. 

shown*4,5’  to  possess  nonspecular  features  at  critical 
Lamb  mode  angles. 

It  was  realized  that  the  distinct  features  associ¬ 
ated  with  nonspecular  reflection  could  be  used  to 
characterize  the  elastic  properties,  or  any  derivation 
therefrom,  of  solid  materials.  These  features  can  con¬ 
stitute  a  readily  recognizable  signature  that  is  uniquely 
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related  to  small  variations  of  certain  elastic  properties 
of  the  medium  to  be  examined.  Mayer  and  Ngoc16-7’ 
have  demonstrated  the  utility  of  nonspecular  features 
as  a  signature  to  detect  local  fluctuations  of  elastic 
properties.  Their  experimental  results  (see  Fig.  2)  and 
preliminary  calculations  showed  that  the  profile  of  a 
2  MHz  beam  reflected  from  a  water-brass  interface 
underwent  a  substantial  shift  in  intensity  between 
two  sections  of  a  beam  as  the  beam  was  scanned 
across  the  known  location  of  an  elastic  defect. 
Chimenti  and  Adler'8*  have  used  Rayleigh-angle  inci¬ 
dence  for  the  detection  of  small  surface-connected 
fatigue  cracks  in  a  titanium  alloy,  demonstrating  that 
the  presence  of  a  fatigue  crack  leads  to  a  reduction  in 
intensity  over  one  section  of  the  nonspecularly  re¬ 
flected  beam. 

This  paper  is  a  theoretical  investigation  of  the 
relationships  between  characteristics  of  elastic  defects 
and  nonspecular  features  of  a  bounded  beam  re- 


Fi*.  2.  Schlieren  photographs  showing  profile  changes  due  to  an 
elastic  defect  for  a  2  MHz  beam  incident  onto  a  water-brass 
interface  at  the  Rayleigh  angle,  (a)  no  defect;  (b)  with  defect. 
(After  ref.  6.) 
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fleeted  from  a  liquid-solid  interface  at  the  Rayleigh 
angle.  It  is  then  intended  to  serve  as  a  theoretical 
basis  for  interpretation  of  Rayleigh-angle  nonspecu¬ 
larly  reflected  profiles,  when  they  are  used  as  a 
signature  to  characterize  a  localized  surface  elastic 
defect. 


2.  NONSPECULAR  REFLECTION  AT 
RAYLEIGH  ANGLE 

Mathematically,  a  bounded  ultrasonic  beam  can 
be  described  as  a  superposition  of  an  infinite  number 
of  plane  waves  of  different  amplitude  which  have  the 
same  frequency  but  are  incident  at  slightly  different 
angles  centered  about  the  principal  incident  direc¬ 
tion.  In  Fig.  1,  an  ultrasonic  beam  bounded  in  the 
(x.  z)  plane  and  uniform  in  the  y  dimension  is  shown 
to  be  incident  at  an  angle  9t,  which  is  assumed  to  be 
the  principal  incident  direction.  The  beam  is  char¬ 
acterized  by  angular  frequency  u  =  2irf  and  con¬ 
tained  within  an  effective  width  of  2w,  which  is 
projected  onto  the  interface  along  the  x  axis  as  2w’0. 

It  was  shown  through  Fourier  analysis'8*  that  the 
sound  field  of  an  incident  beam  can  be  uniquely 
determined  at  any  point  in  the  (x,  z)  plane  if  its  field 
distribution  is  given  in  any  plane.  Consequently,  once 
the  sound  field  is  determined  in  the  plane  ;  =  0.  a 
well-defined  incident  beam  can  be  specified  anywhere 
by  the  following  Fourier  integral  transform  pair: 

Uinc(x.z)  =  (2*)-'fk‘  +  n/WoV(kx) 

Jk,  -  lt/w0 

■exp[i(xkx  +  zk.)]dkx  (1) 

v(  kx )  =  /  ^inc ( x , 0 )exp(  -  ixkx )  dx  (2) 

In  (1)  and  (2),  k  is  the  wave  vector  variable 
whose  x  and  z  components  are  kx  and  k.\  the  wave 
vector  k,  =  (w/t>)  sint?  describes  the  central  direc¬ 
tion  of  propagation  of  the  incident  beam;  V(kx)  is 
the  amplitude  of  the  constituent  plane  waves; 
l/,nc(x,0)  is  the  field  distribution  of  the  incident 
beam  in  the  plane  z  =  0;  and  the  time-dependent 
factor,  exp(-/wt),  is  suppressed  for  convenience.  In 
all  calculations  described  in  this  paper,  the  incident 
beam  is  assumed  to  have  a  Gaussian  distribution, 
which  is  characterized  at  the  plane  z  =  0  by 

t/inc(-x>0)  =  cxp[-(x/wo)2+ixAx] 


(3) 
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Hence,  from  (2),  the  specific  form  of  V(kx)  is 

V(kx)  -  ir'/2w0exp[-(/cJ[  -  fc,)2( h>0/2)2]  (4) 

The  mathematical  representation  that  describes 
the  incident  beam  by  (1)  can  be  extended  to  represent 
a  bounded  ultrasonic  beam  reflected  from  a  liquid- 
solid  interface.  In  the  same  manner,  the  reflected 
beam  is  a  superposition  of  constituent  plane  waves, 
which  are  individually  reflected  according  to  the  plane 
wave  reflection  coefficient  R(kx)  associated  with  the 
liquid-solid  structure.  In  other  words,  the  sound  field 
of  the  reflected  beam  can  be  represented  by 

U(x^)^(2irr'fkl  +  V/^R(kx)V(kx) 

Jk,-ir/w0 

exp  [i(xkx  +  zk,)\dkx  (5) 

where  R(kx)  is  given<2)  by 

r(k  )  (kf-2klf  +  4k2xK,Kd-pk*Kd/K 

( k]  - 2k\ )'  +  4 k\K,Kd  +  pk*Kd/K 

in  which 

p-p'/p" 

k  —  (u/v)(\  +  ia/2ir) 
kd  =  (u/vd){\  +  iad/2ir) 

)( 1  +  ia,/2it ) 

K.-(kl-kl)vl 

K,-(k)-kl)'n  (7) 

The  quantities  a ,  ad ,  and  a,  in  (7)  denote  attenuation 
per  wavelength  for  sound  waves  in  liquid  and  for 
longitudinal  and  shear  waves  in  solid,  respectively. 
The  integral  in  (5)  can  be  carried  out  to  provide  the 
distribution  profile  of  the  sound  field  of  the  reflected 
beam  for  any  incident  angle.  For  example,  the  inten¬ 
sity  profile  of  a  2  MHz,  20  mm  Gaussian  beam 
reflected  from  a  water-stainless  steel  interface  at  the 
Rayleigh  angle  is  presented  in  Fig.  3,  where  the 
profile  (solid  curve)  is  seen  to  have  primarily  two 
sections,  with  the  higher  one  being  laterally  displaced 
to  the  right  and  the  lower  one  to  the  left  of  the  center 


*H» 


Fig.  3.  Variations  of  reflected  profile  as  t>x  is  changed  with  vL  and 
density  is  kept  constant.  (p"~  7.94,  vL  —  5.84  km/s,  and  S,  — 
corresponding  to  cs  -  3.13  lun/s). 

of  the  beam.  This  is  a  typical  nonspecularly  reflected 
profile  having  two  sections  separated  by  an  intensity 
minimum  and  a  trailing  field  away  from  the  insoni- 
fied  area. 


3.  APPLICATION  OF  RAYLEIGH- ANGLE 
NONSPECULAR  REFLECTIVITY  TO 
DEFECT  CHARACTERIZATION 

The  very  distinct  features  produced  by  the  non- 
specular  effects  present  themselves  as  an  effective 
means  to  characterize  elastic  properties  associated 
with  the  material  medium  under  investigation,  pro¬ 
vided  small  fluctuations  in  the  elastic  properties  can 
be  related  to  discernible  changes  in  the  nonspecular 
features  of  the  reflected  profile.  In  order  to  determine 
whether  these  changes  are  discernible,  or  if  any  exist 
at  all,  reflected  profiles  were  calculated  for  Rayleigh- 
angle  incidence  as  the  key  elastic  properties  of  stain¬ 
less  steel  were  varied.  Starting  with  the  initial  selec¬ 
tion  of  elastic  parameters,  p  -  7.94,  vL  -  5.84  km/s, 
and  -  3.13  km/s,  each  of  these  parameters  was 
varied  while  keeping  the  others  fixed.  It  was  found 
that  the  nonspecular  features  were  quite  insensitive  to 
variation  of  p  and  vL.  However,  even  small  fluctua¬ 
tions  in  vs  resulted  in  substantial  changes  in  the 
nonspecularly  reflected  profile. 

Figure  3  shows  the  results  of  calculations  where 
the  value  of  the  shear  wave  velocity  was  changed 
from  the  assumed  value  of  vs  -  3.13  km/s  to  3.09 
and  3.05  km/s,  respectively.  It  is  seen  that  a  change 
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Fig.  4.  Percentage  change  tn  as  a  function  of  percentage  change 
in  v,  and  vs.  respectively. 

of  about  1.5%  in  the  value  of  vs  can  produce  substan¬ 
tial  changes  in  the  reflected  beam  profile  as  the 
higher  section  shifts  to  the  left  while  the  intensity 
minimum  becomes  less  pronounced.  These  results 
indicate  that  the  shear  wave  velocity  of  the  solid 
medium  has  a  dominant  influence  in  determining  the 
nonspecular  features.  This  can  be  seen  from  Fig.  4, 
where  the  percentage  changes  in  vR  are  plotted  as 
functions  of  percentage  changes  in  v,  and  es,  demon¬ 
strating  that  vR  is  strongly  related  to  vs  but  is  almost 
independent  of  small  changes  in  vL. 

It  can  therefore  be  established  from  the  preced¬ 
ing  results  that,  in  applying  Rayleigh-angle  nonspecu¬ 
lar  reflectivity  to  elastic  defect  characterization, 
defects  can  be  reasonably  specified  by  the  value  of 
the  shear  velocity  only.  The  simplest  type  of  elastic 
defects  can  then  be  specified  by  setting,  for  example, 

_  1 3.23  km/s  -  w0  «  x  «  0 
5  \  3.13  km/s  otherwise 

This  stepwise  defect,  which  is  of  half-beamwidth  size 
and  centered  at  x  -  -  tv0/2,  is  schematically  repre¬ 
sented  by  the  curve  of  a  step  function  below  the  x 
axis  in  Fig.  5.  As  the  beam  impinges  upon  the  surface 
area  which  contains  such  a  defect,  the  nonspecularly 
reflected  beam  is  expected  to  have  a  quite  different 
profile.  Figure  5  shows  a  sharp  increase  in  intensity 
in  the  left  section  of  the  reflected  beam,  where  the 
elastic  defect  is  present.  With  -  3.23  km/s,  vs 
changes  by  3%  as  the  peak  intensity  of  the  left  section 
of  the  beam  becomes  four  times  stronger.  However, 


Fig.  5.  Changes  in  reflected  profile  (upper  curves)  due  to  three 
stepwise  changes  in  t'x  (lower  curves). 

with  a  much  smaller  percentage  change,  1  %  ( vs  =  3. 1 7 
km/s),  the  change  in  the  peak  intensity  of  the  left 
section  is  readily  observable. 

With  elastic  defects  simply  described  by  (8),  the 
above  results  can  serve  to  demonstrate  the  possibility 
of  applying  nonspecular  effects  of  Rayleigh-angle 
reflection  to  defect  characterization.  Another  specifi¬ 
cation  of  elastic  defects,  which  is  probably  more 
realistic  than  that  given  by  (8),  can  be  achieved  by 


Ftg.  6.  Changes  in  reflected  profile  due  to  a  smoothly  varying 
elastic  defect  (upper  curves).  Variations  of  t's  are  assumed  to  be 
sinusoidal  (lower  curves).  Profile  change  shown  in  Fig.  5  is  re¬ 
drawn  for  comparison. 
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representing  the  shear  velocity  vs  in  the  form 

1 3.18  +  0.05  sin[2ff(x  +  xc)/s  +  it/ 2] 

Vs  =  /  for  xc  -  s/2  <  x  «  xc  +  s/2 

\3.13  for  x  >  xc  +  s/2  and  x  <  xc  -  s/2 

(9) 

An  elastic  defect,  as  modeled  by  (9),  has  size  s,  and 
its  center  is  located  at  x  =  xc.  A  defect  of  type  (9)  is 
plotted  in  Fig.  6  along  the  x  axis,  where  s  =  h>0  and 
xc  =  -  *t'0/2.  With  this  defect  present.  Fig.  6  shows 
the  profile  of  a  2  MHz,  20  mm  beam  as  it  is  reflected 
at  the  Rayleigh  angle  from  the  water-stainless  steel 
interface.  The  left  section  of  the  beam  is  shown  to  be 
affected  by  the  presence  of  the  defect,  having  a  much 
higher  intensity.  The  resulting  reflected  profile  now 
does  not  show  any  abrupt  variation  as  shown  in  Fig. 
5  when  a  stepwise  defect  is  considered. 


4.  DETERMINATION  OF  DEFECT  SIZE  AND 
LOCATION 

The  above  calculations  have  illustrated  how  non- 
specular  reflectivity  at  the  Raleigh  angle  can  be 
utilized  to  detect  an  elastic  defect.  In  other  words,  the 
nonspecularly  reflected  profiles  are  used  as  signatures 
that  can  determine  whether  such  a  defect  exists.  This 
application  can  be  made  more  effective  if  informa¬ 
tion  derivable  from  these  signatures  can  point  out 
where  a  defect  is  located,  how  large  it  is,  and  to  what 
extent  the  elastic  property  deviates  from  the  standard 
value.  The  following  calculations  are  intended  to 
illustrate  how  nonspecularly  reflected  profiles,  when 
used  as  defect  signatures,  are  sensitive  to  variations  in 
defect  size,  location,  and  elastic  deviation.  In  these 
calculations,  the  elastic  defect  under  consideration  is 
the  smoothly  varying  type  described  by  (9). 

First,  three  reflected  profiles  were  calculated  for 
three  locations  of  the  elastic  defect  with  its  size  being 
kept  the  same,  s  =  w0.  The  results  are  shown  in  Fig.  7 
for  xc  =  —  w0/2.  0,  and  w0/2.  In  Fig.  7,  as  the  defect 
is  moved  from  left  to  right,  there  appears  an  addi¬ 
tional  intensity  peak  that  moves  with  it.  Since  the 
elastic  defect  causes  the  elastic  property  to  deviate 
from  that  corresponding  to  Rayleigh-angle  incidence, 
it  is  apparent  that  the  part  of  the  incident  beam 
impinging  upon  the  defect  area  undergoes  a  specular 


Ft*.  T  Variation  of  defect  characterizing  profile  as  a  function  of 
defect  location  with  defect  centered  (a)  at  x/w0-  -0.5;  (b)  at 
v/h0  -  0;  (c)  at  x/wB-  +0.5.  Defect  size  and  deviations  as  in 
Fig.  6. 


reflection.  This  specularly  reflected  beam  section  con¬ 
firms  the  existence  of  the  defect  and  also  identifies 
where  the  defect  is.  In  addition,  the  calculation  re¬ 
sults  presented  in  Fig.  5  lead  to  the  observation  that 
the  more  the  elastic  defect  deviates  from  the  standard 
value,  the  stronger  the  intensity  of  the  specular  sec¬ 
tion  is. 


Ftg*.  Variation  of  defect  characterizing  profile  as  a  function  of 
defect  size.  Variation  in  defect  size  (lower  curve)  causes  changes  in 
beam  profile  (upper  curve). 
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Next,  the  reflected  profiles  were  calculated  for 
different  sizes  of  the  elastic  defect,  s  -  *v0,  w0/2,  and 
w0/4.  The  location  of  the  defect  is  determined  by 
«  0  and  remains  the  same  for  all  three  calcula¬ 
tions.  The  results  are  presented  in  Fig.  8  where  the 
defect  size  is  characterized  along  the  position  axis.  It 
is  established  in  this  figure  that  the  smaller  the  size  of 
the  defect  becomes,  the  smaller  the  width  of  the 
specular  beam  section  is.  In  other  words,  the  size  of 
the  defect  can  be  estimated  by  the  width  of  the 
specular  beam  section. 


5.  CONCLUSIONS 

In  summary,  the  above  calculation  results  have 
demonstrated  that  nonspecular  reflectivity  at  the 
Rayleigh  angle  can  be  used  to  characterize  a  localized 
surface  elastic  defect.  In  this  application,  the  non- 
specularly  reflected  profile  serves  as  a  signature  of 
the  elastic  defect,  whose  size,  location,  and  deviation 
of  the  shear  velocity  can  be  estimated  by  the  width, 
position  and  intensity  level  of  the  specular  beam 
section  associated  with  the  defect. 
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Measurements  of  harmonic  profiles 
of  a  bounded  ultrasonic  beam  in  a 
liquid  medium 

D.D,  MCLENNAN,  T.D.K.  NGOC  and  W.G.  MAYER 


The  acoustic  field  of  a  bounded  finite-amplitude  beam  is  governed  by  competing 
mechanisms  such  as  medium  non-linearity,  interference,  and  attenuation.  An  analytic 
determination  of  the  distribution  of  sonic  energy  in  the  individual  harmonic 
components  is  often  not  reliable;  consequently  an  experimental  mapping  of  those 
harmonic  components  may  be  more  appropriate.  An  experimental  method  is  presented 
here  which  involves  a  miniature  hydrophone  probe  which  has  been  frequency  calibrated 
for  each  harmonic  component  using  a  light  diffraction  technique.  Measurements  (made 
across  the  sound  field)  of  the  acoustic  amplitude  of  the  first  three  harmonic  components 
are  presented. 

KEYWORDS:  ultrasonics,  transducer  calibration,  hydrophones 


Introduction 

Acoustic  wave  propagation  is  an  inherently  non-linear 
process.  As  a  consequence  of  this  non-linearity ,  an 
originally  sinusoidal  waveform  distorts  toward  a  sawtooth 
shape  as  it  travels  through  a  medium.  Examination  of  the 
spectral  content  of  the  distorted  waveform  reveals  that  this 
distortion  corresponds  to  the  self-generation  of  harmonics 
within  the  sound  beam.  Thus,  the  harmonic  content 
represents  a  quantitative  measure  of  the  waveform 
distortion. 

Many  ultrasonic  experiments  are  conducted  under  con¬ 
ditions  conducive  to  the  production  of  significant  wave 
distortion  (that  is,  some  combination  of  high  initial  intensity 
and  long  pathlengths).  Moreover,  numerous  evaluation  tech¬ 
niques  using  ultrasonics  are  based  on  the  assumption  that 
the  probing  beam  contains  only  one  frequency.  If,  however, 
sonic  distortions  are  present,  beam  reflectivity,  absorption, 
intensity  distribution  within  the  beam  and  other  parameters 
may  change  sufficiently  so  that  the  evaluation  of  the 
received  beam  may  lead  to  erroneous  results  if  one  does  not 
take  account  of  the  multiple  frequency  components  present 
in  the  originally  sinusoidal  beam. 

Rather  than  attempting  to  determine  analytically  the  distri¬ 
bution  of  ultrasonic  energy  as  well  as  the  degree  of  distortion 
present  in  all  areas  of  the  sound  field,  an  often  easier  and 
more  reliable  method  is  to  probe  this  sound  field  with  a 
measuring  device  sensitive  to  the  harmonic  constituents. 

Measurements  of  the  harmonic  distribution  in  an  acoustic 
beam  have  been  made1  •*  showing  the  profiles  of  the  funda¬ 
mental  and  some  harmonics  at  several  propagation  distances. 


The  authors  ara  In  tha  Physics  Department,  Georgetown  University. 
WaWington  DC  20067,  USA.  Paper  received  16  January  1962. 
Revised  11  June  1962. 


But,  to  obtain  quantitative  comparisons  between  the 
frequency-content  distribution  of  the  sound  fields  at  dif¬ 
ferent  locations,  it  is  necessary  to  make  measurements  with 
a  probe  whose  frequency  response  is  accurately  calibrated. 

In  this  paper,  experimental  investigation  of  non-linear 
harmonic  generation  involves  quantitative  measurements  of 
the  sound  field  of  a  low-MHz  ultrasonic  beam  propagating 
in  water  using  a  calibrated  miniature  hydrophone. 

Measurements  were  made  in  a  region  of  the  sound  field  that 
was  useful  experimentally,  as  well  as  in  an  area  where  all 
competing  effects  such  as  non-linearity,  attenuation,  and 
interference  are  significant.  In  the  following,  the  experimen¬ 
tal  procedure  is  described  in  detail  and  the  measured 
profiles  of  the  fust  three  harmonics  are  presented. 

Experimental  procedure 

The  experimental  procedure  is  centred  around  two 
techniques: 

( 1 )  Mapping  of  the  sound  field  by  an  electro-mechanical 
probe. 

(2)  Acousto-optic  interaction  for  calibrating  the  probe. 

The  experimental  set-up  is  therefore  designed  to  accommo¬ 
date  both  techniques.  Because  the  calibration  procedure 
is  rather  involved,  it  will  be  described  separately. 

Experimental  methodology 

Of  the  known  methods  currently  in  use,5  the  two  techniques 
most  often  mentioned  for  mapping  bounded  acoustic  beams 
are  miniature  hydrophones  and  diffraction  of  light  by 
sound.  The  miniature  hydrophone  is  suited  for  probing  a 
small  area,  but  it  does  not  have  a  linear  frequency  response.4 
On  the  other  hand,  the  light  intensity  of  a  diffraction  order 
is  independent  of  frequency  but  represents  an  integration 
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of  the  acoustic  field  across  the  sound  beam.  The  experimen¬ 
tal  technique  used  here,  employing  a  miniature  hydrophone 
to  map  the  sound  field  and  a  light  diffraction  set-up  to 
calibrate  the  hydrophone,  makes  use  of  the  major  strengths 
of  these  two  techniques. 

The  specific  experimental  set-up  is  illustrated  in  Fig.  1 .  The 
sound  beam  is  generated  along  the  z-axis  and  the  translation 
of  the  probe  is  along  the  x-axis,  across  the  sound  beam.  The 
electrical  signal  from  the  probe,  representing  the  local 
sound  field,  is  connected  to  a  spectral  analyser  to  select  the 
desired  harmonic  component,  the  sound  absorber  is 
installed  behind  the  probe  to  ensure  that  the  ultrasonic 
beam  consists  of  progressive  waves  only.  The  light  diffrac¬ 
tion  system  is  also  presented  in  Fig.  1 .  The  combination  of 
light  source,  lenses,  and  aperture  produces  a  collimated 
light  beam  passing  perpendicularly  to  the  sound  beam  and 
forms  a  diffraction  pattern  on  the  face  of  the  photodiode. 
The  lenses  and  aperture  were  adjusted  so  as  to  maximize  the 
spatial  separation  of  the  diffraction  orders. 

The  source  of  the  ultrasonic  field  was  an  air-backed  quartz 
transducer  driven  at  3  MHz  and  having  an  effective  radiating 
diameter  of  2.0  cm.  The  probe  used  was  a  miniature 
hydrophone  manufactured  by  Mediscan  Incorporated  with 
a  circular  detecting  area  1 .0  mm  in  diameter  embedded  in  a 
hypodermic  needle  1 .3  mm  in  diameter.  The  detecting  area 
of  the  probe  is  thus  comparable  to  the  sound  wavelength 
and  would  therefore  represent  only  a  minimal  disturbance 
to  the  sound  field.  In  addition,  the  probe  was  inserted 
through  a  piece  of  sound  absorbing  rubber  positioned  in 
front  of  the  electrical  connection  of  the  probe,  thus 
eliminating  the  possibility  of  establishing  a  significant 
standing-wave  pattern  between  transducer  and  probe. 

Calibration 

The  measuring  probe  was  calibrated  using  the  diffraction 
of  light  by  sound.  The  interaction  of  light  with  high- 
frequency  sound  manifests  itself  in  the  form  of  a  diffraction 
pattern  with  the  light  intensity  of  the  nth  order  being  given 
by5 


In  =  0) 

where  J„(v)  is  the  nth-ordered  Bessel  function  of  the  first 
kind.  The  argument  v,  the  Raman-Nath  parameter,  of  the 
above  Bessel  function  is  defined  as 

v  *=  kp0L  (dM/dp)s  (2) 


FH.1  Schematic  of  the  experimental  wt-up.  The  sound  beam 
trevele  in  the  /  direction  and  the  light  beam  in  the  *  direction 


Fig.  2  Plot  of  Raman-Nath  parameter  against  probe  output  for 
3,  6.  and  B  MHz.  The  parameter  v  is  dimensionless  and  the  sound 
amplitude  is  expressed  in  arbitrary  units 


where  k  is  the  wavenumber  of  light,  p0  the  maximum 
acoustic  pressure,  L  the  interaction  length,  and  (dju/dp)s 
the  partial  derivative  of  the  refractive  index  with  respect 
to  pressure  in  adiabatic  conditions.  This  partial  derivative 
is  identified  as  the  piezo-optic  coefficient.  When  a  bounded 
sound  beam  is  considered,  the  value  of  p0  in  (2)  must  be 
replaced  by  the  maximum  acoustic  pressure  averaged  over 
the  interaction  length. 

Within  the  narrow  frequency  range  under  consideration, 
the  piezo-optic  coefficient,  and  hence  the  Raman-Nath 
parameter,  is  independent  of  sound  frequency.6  As  a  result, 
(2)  implies  that  there  exists  a  linear  relationship  between 
the  Raman-Nath  parameter  v  and  the  acoustic  pressure  p0 , 
which  is  a  measure  of  sound  amplitude.  In  turn,  v  can  be 
determined  from  the  measured  light  intensity  using  (1). 

The  output  of  the  measuring  probe,  which  is  considered  to 
consist  of  the  miniature  hydrophone  and  the  spectrum 
analyser,  is  calibrated  against  v  derived  from  the  known 
behaviour  of  /0J(v)  by  measuring  the  percentage  light 
intensity  of  the  zeroth  order.  With  reference  to  Fig.  1 ,  the 
calibration  was  carried  out  with  the  light  beam  intercepting 
the  sound  field  between  the  transducer  and  the  hydrophone, 
which  was  positioned  about  3  mm  from  the  transducer  face. 
In  this  manner,  the  diffracted  light  intensity  and  the  probe 
output  were  simultaneously  measured  at  different 
transducer  driving  voltages  for  3, 6,  and  9  MHz. 

The  values  of  v  were  derived  following  the  above  procedure 
and  were  then  plotted  against  the  probe  output  u  illustrated 
in  Fig.  2.  It  can  be  seen  that  the  relationship  between  v  and 
the  probe  output  is  linear  for  all  frequencies.  The  dopes  of 
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the  three  lines  in  Fig.  2,  then,  represent  the  specific  response 
of  the  measuring  probe  to  each  frequency.  Relative  to  the 
fundamental  frequency,  this  response  can  be  quantitatively 
described  in  terms  of  the  ratios  1.0:2.9:3.5  for  3, 6,  and 
9  MHz,  respectively. 

Measurement  of  harmonic  profiles 

The  above  procedure  was  used  to  measure  the  local  sound 
field  as  the  sound  waves  travelled  in  water.  Throughout  the 
experiment  the  launched  intensity  was  kept  the  same  by 
setting  the  sound  amplitude  at  the  beam  centre  to 
correspond  to  v  =  2.4,  which  is  the  first  zero  of  Jo(v), 
describing  the  light  intensity  of  the  zeroth  order.  This  was 
also  used  as  a  procedure  to  check  whether  standing  waves 
existed  in  the  sound  field,  because  the  light  intensity  of 
the  zeroth  order7  would  not  behave  like  /02(v)  and  would 
not  go  to  zero  at  v  -  2.4  in  the  presence  of  standing  waves. 
In  this  experiment,  the  absence  of  significant  standing  waves 
was  continually  confirmed  by  this  procedure. 

Measurements  were  made  across  the  sound  beam  at  six 
different  probe-to-source  distances  for  3, 6,  and  9  MHz. 

Fig.  1  indicates  that  the  position  of  the  probe  across  the 
beam  is  described  by  the  x-axis  and  the  probe-to-source 
distance  by  the  z-axis.  As  a  result,  three  sets  of  beam  pro¬ 
files  were  obtained  describing  the  variation  of  beam 
profiles  as  a  function  of  probe-to-source  distance  under  the 
influence  of  competing  mechanisms:  non-linearity,  attenua¬ 
tion,  and  interference.  These  are  presented  in  Figs  3-6. 

Figure  3  shows  typical  beam  profiles  of  three  harmonic 
components  measured  at  the  same  distance,  z  =  21  cm.  The 
vertical  scale,  in  arbitrary  units,  represents  the  measuied 
sound  amplitude,  which  has  been  adjusted  in  accordance 
with  the  calibration  results  obtained  earlier.  The  lines  con¬ 
necting  the  points  are  meant  to  be  a  visual  aid  and  do  not 
represent  a  fit  to  the  data. 


F*.  3  Bvprw ntetiw  profile*  of  ttw  first  thrw  harmonics  at 
z  -  21  cm 


Fig.  4  Variation  of  beam  profile  along  the  i  direction  for  the 
fundamental  frequency  component 


*  [cm] 


Fig.  6  Variation  of  beam  profile  along  the  t  direction  for  the 
6  MHz  component 


Figures  4-6  present  the  variation  of  beam  profiles  of  3, 6, 
and  9  MHz  components  along  the  z-axis.  The  values  of 
probe-to-source  distance  range  from  15cmto25cmin2cm 
increments.  These  probe-to-source  distances  were  chosen 
for  lnvestiption  in  order  to  have  measurements  made  in  a 
range  where  effects  due  to  non-linearity,  attenuation,  and 
interference  were  believed  to  be  significant  while  all  har¬ 
monic  components  of  interest  were  still  measurable. 

Some  salient  features  of  the  beam  profiles  can  be  observed 
in  Figs  4-6.  The  oscillatory  variation  in  the  maximum 
amplitude  of  the  3  MHz  profiles  shown  in  Fig.  4  points 
to  the  dominance  of  the  interference  process  in  determining 
the  general  profile  shape.  In  Figs  S  and  6,  the  presence  of 
the  medium  non-linearity  mechanism  is  dearly  seen  in  the 
steady  growth  of  the  6  and  9  MHz  profiles.  Particularly  in 
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Fig.  5,  at  z  ■  23  and  25  cm,  the  growth  of  the  6  MHz  beam 
profiles  seems  to  be  levelling  off,  indicating  a  possible  balance 
between  the  damping  effect  of  attenuation  and  the  growth 
of  the  harmonic  due  to  medium  non-linearity.  For  the 
near -Gaussian  initial  profile  produced  by  the  3  MHz  trans¬ 
ducer  under  investigation,  it  is  observed,  in  general,  that  the 
profiles  of  all  harmonic  components  do  not  show  any 
appreciable  beam -spreading  effect. 

Conclusion 

The  above  results  establish  that  the  employed  experimental 
method,  which  combines  miniature  hydrophone  measure¬ 
ment  and  a  light  diffraction  technique,  allows  for  a  quanti¬ 
tative  comparison  of  different  harmonic  components  of  a 
finite-amplitude  bounded  acoustic  beam  travelling  in  a 
liquid  medium.  From  these  results  one  can  draw  conclusions 
about  the  actual  frequency-content  of  an  ultrasonic  beam 
which  has  travelled  in  a  liquid  medium  before  it  is  used  to 
evaluate  the  properties  of  a  sample  to  be  tested.  This  know¬ 
ledge  appears  to  be  important  in  all  cases  where  the 
evaluation  technique  is  frequency  dependent. 
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